Recursive graphs with small-world scale-free properties 
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We discuss a category of graphs, recursive clique trees, which have small-world and scale-free 
properties and allow a fine tuning of the clustering and the power-law exponent of their discrete 
degree distribution. We determine relevant characteristics of those graphs: the diameter, degree 
distribution, and clustering parameter. The graphs have also an interesting recursive property, and 
generalize recent constructions with fixed degree distributions. 
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In a complex system a large number of agents interact 
showing a cooperative behavior which does not depend 
only on the individual features of its parts, but also on 
its structure and on the different sort of relations which 
can be established between them and the environment. 
This global behavior allows the system to attain certain 
achievements without the presence of an administrative 
hierarchy or a central control mechanism. A swarm of 
bees, an ant colony, companies which supply a big city 
(water, electricity, telephone, etc.), even biological mech- 
anisms or social relationships are all complex systems 
where global patterns emerge from the interaction of a 
large number of similar elements. In recent years, there 
has been an increase in understanding of such complex 
systems in terms of networks, modeled by graphs com- 
posed of vertices and edges, where vertices represent the 
basic elements and edges their interactions. The World 
Wide Web (WWW), Internet, transportation systems, 
and many biological and social systems have been char- 
acterized by small- world scale-free networks 0, Q] which 
have a strong local clustering (nodes have many mu- 
tual neighbors), a small diameter (maximum distance be- 
tween any two nodes) and a distribution of the number of 
nearest neighbors for each node according to a power law 
(scale- free property), see, for example, Refs. HSHH. 
To model these networks and their growth very often 
stochastic models and methods from statistical physics 
have been considered; see Ref. Q. However, the use of 
exact deterministic models allows a quick determination 
of the relevant parameters of the associated graph that 
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may be compared with experimental data from real and 
simulated networks. Previous work, for example, consid- 
ered deterministic small-world networks comparable to 
those obtained stochastically by Watts and Strogatz 
node replacement and node addition methods to produce 
small- world and scale-free networks from a low diameter 
"backbone" network [jj, specific recursive scale- free con- 
structions with fixed degree distribution s IToL ITlL Il3| 
and scale- free trees (without clustering) jlj]. In this pa- 
per we present a recursive graph construction which pro- 
duces scale-free small- world networks with an adjustable 
clustering and such that the parameter of the power law 
associated with the degree distribution, 7, takes values 
between 2 and l + ln3/ln2 = 2.58496. It has been shown 
that the values for the scaling exponent of most real tech- 
nical networks 0] are in this range. Moreover, the ad- 
dition of one-degree nodes to our construction allows a 
model which is not far from the real topology of Inter- 
net at the autonomous system level, as is described in 
Ref. [3. 

Therefore, we introduce here a deterministic exact net- 
work model of complex systems for which we can adjust 
both the clustering parameter and the power-law expo- 
nent, providing a complementary tool to stochastic ap- 
proaches. 

Recursive clique tree K(q,t), q > 2. Definition. — a 
complete graph Kj, (also referred in the literature as q 
clique; see Ref. [lj]) is a graph without loops whose q 
vertices are pairwise adjacent. The recursive clique tree 
K(q,t) is the graph constructed as follows: For t = 0, 
K(q, 0) is the complete graph K q (or q clique). For t > 1, 
K(q,t) is obtained from K(q,t — 1) by adding for each 
of its existing subgraphs isomorphic to a g-clique a new 
vertex and joining it to all the vertices of this subgraph 
JiirtS^ fflfor the case q — 3). Then, at t = 1, K(q, 1) 
results in the complete graph with q+l vertices, K q+ \ 1 
and at t — 2 we add q + l new vertices each of them 
connected to all the vertices of one of the q cliques K q 
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FIG. 1: (Color online) First stages of a growing recursive 
clique tree K(3,t). K(3,0) is the complete graph or clique 
K3 (the triangle with thick edges). At each step a vertex is 
connected to each of existing K3. Therefore at t = 1 only 
one vertex is added, but at t — 2, four new vertices should be 
connected to the four different cliques K3. 




FIG. 2: (Color online) Recursive construction of K(2,i). We 
glue a native (or initial) clique of a K(2,t—1) on each 2-clique 
of Kz to construct K(2,t). In this figure we obtain K(2, 3) 
from three copies of K(2, 2) as follows: The graph at the left is 
K(2, 2) with its native 2-clique (the thick edge). Three copies 
of K(2, 2) are glued by their natives 2-cliques to each of the 
three 2-cliques of K3 (center) resulting in K(2,3) (right). 



(subgraphs of K q+ i), and so on. 

This algorithm produces a complex growing graph with 
a tunable parameter q which controls all its relevant char- 
acteristics. In the particular case q = 2 we obtain the 
same graph as in Ref. fill ]. However, our family is infi- 
nite as q can take any natural value starting from 2. 

Notice that although we call the graph a recursive 
clique tree, the graph contains numerous cycles and hence 
is not a tree in the strict sense. Recursive clique-trees 
K(q, t) are a natural generalization of trees (if one con- 
siders the case q = 1, then we obtain a subfamily of the 
"usual" trees). Similar constructions in which new ver- 
tices are joined to every vertex of a given q clique have 
been considered, in another context, in Ref. |l8| and 
termed u k trees." 

Recursive construction of K(q, t). There exists another 
interesting way to construct K (q, t), which clearly shows 
the recursive structure of such networks. We call native 
clique of K(q, t) its initial q clique at t = 0. Then K(q, t) 
is constructed as follows: 

At step 0, we have the native clique K q . For any step 
t > 1, K(q,t) is constructed as follows: Consider a (q + 
1) clique, then every subgraph of it isomorphic to a q 
clique (there are q + 1 such cliques) is a native clique of 
a K(q,t- 1); see Fig.0 

Size and order of K(q,t). Table [I] gives the number 



TABLE I: Number of new edges added to the recursive clique 
tree K(q,t) at each step t and total number of complete 
graphs K q at this step. 

Step (t) New edges Number of K q 



1 q q + 1 

2 q(q + l) q(q + l) + (q + l) = (q+l) 2 

3 q(q + lf q(q+l) 2 + (q + l) 2 = (q + 1) 3 
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TABLE II: 


Distribution of vertices and degrees 


for the recur- 


sive clique tree K(q, t) at each step t. 
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of new edges added to the tree at each step and the to- 
tal number of K q at this step. Therefore we can easily 
compute the total of edges at step t, 

\E\ t = ^+q£(q+iy= q -^ + (q+iy-l 

i=0 

(1) 

The distribution of vertices and degrees is given in Ta- 
bled! 

The maximum degree at step i is 
q — <? q % — 1 

Aj = + q = ~ + q ~ 1 2 

1 - q q-1 

Adding up the number of vertices gives the result 

^=^(g + l) J + (g + l)= {<1 + 1)4-1 +q (3) 

H q 

The average degree is then 

j _ 2\E\t _ q[q 2 ^q + 2(g+l) t -2} 

* N t ( g+ l)*_l + g2 {! 
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[for q = 2, it is 4/(1 + 3 1 -*)] 

Degree distribution. The degree spectrum of the graph 
is discrete: at time t, the number N(k, t) of vertices of 

degree k = q, 2q, q 2 + 2q, q 3 +q 2 +2q, . .., g t_2 +g t_3 H h 

g 2 +2g,g' _1 +g t ~ 2 H \-q 2 + 2q is equal to (g+l)*" 1 , (q+ 

l)*- 2 , (g+l)*- 3 , (<z+l)*- 4 , . . • , (9+1), (?+l), respectively. 
Other values of the degree are absent. Clearly, for the 
large network, N(k,t) decreases as a power of k, so the 
network can be called "scale free." Spaces between de- 
grees of the spectrum grow with increasing k. Therefore, 
to relate the exponent of this discrete degree distribution 
to the standard 7 exponent of a continuous degree distri- 
bution for random scale free networks, we use a cumula- 
tive distribution P cum (k) = Y,k'>k N(k',t)/N t ~ A: 1-7 . 

Here k and k! are points of the discrete degree spec- 
trum. For a degree 
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there are (q + 1)' _1 vertices with this exact degree. 

We count now how many vertices have this and a 
higher degree. From the distribution 



z-i 



£i\T(fc',t) = 5>+l) p + (? + !) = 
P =i 
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As the total number of vertices at step t is given in Eq. © 
we have 



q( 



t-i 



1 



q-i 



l) 



1 -_^ + <z_(g + i)<-i + , 



(g+D'-l + q ( q+ l)t-l + q 2 

Therefore, for t large 

= (Q + 1)'-* 



and 



7 w 1 + 



ln(q+ 1) 
In q 



(5) 



so that 2 < 7 < 2.58496. 

Also, notice that when t gets large, the maximal degree 



of a vertex is roughly q* _1 ~ j\j^ n i/ ln ('}+ 1 ) _ ^ 

Clustering distribution. The clustering coefficient C (x) 
of a vertex x is the ratio of the total number of exist- 
ing connections between all k its nearest neighbors and 
k(k— 1)/2, the number of all possible connections between 
them. The clustering of the graph is obtained averaging 
over all its vertices. In what follows Ni(x) will denote the 
total number of neighbors which x has at step i, /Q(x) 
will denote the number of different q cliques which con- 
tain x at step i, and £%(x) will denote the total number 
of edges among the neighbors of vertex x also at step i. 

Next we compute the clustering parameter for a vertex 
x at any step of iteration. 

Step 1. As x belongs to K q+ i, i.e., the recursive 
clique tree K(q,t) at step 1, it is adjacent to q other 
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vertices [Afi(x) = q] which form a K q . Therefore the 
number of edges among all vertices adjacent to x is 
£\(x) = q(q — l)/2. Moreover, x belongs to q different q 
cliques, i.e., Ki(x) = q. 

Step 2. As the number of q cliques in K q+ \ is q + 1, we 
add this number of new vertices to the construction. But 
of these, q will be forming -each- a different (q + 1) clique 
with x. Therefore N~2{x) — 2q and £%{x) = q(q — l)/2 + 
q(q— 1). Now x will belong to /C 2 (x) = Ki(x)+q(q— 1) = 
q 2 different q cliques. 

Step 3. We add now to the construction K,2(x) new 
vertices which will form each a (q + 1) clique with x. 
£ 3 {x) = £ 2 {x) + K 2 (x)(q - 1) = q(q - l)/2 + g(g - 1) + 
<Z 2 (g-l). /C 3 (x) = /C 2 (x)+/C 2 (x)(g-1) = q 2 + <Z 2 ( (7 -l) = 
g 3 and A/" 3 (x) = (q 3 - l)/(q - 1) + 9 - 1. 

5tep We add /C t _i(x) new vertices forming each a 
(q + 1) clique with x. <£t(x) = £t_i(x) +K.t-i(x)(q— 1) = 

9(9 - l)/2+(«-l) E^l 9 l = (g-l)b/2+q* - 1/(5 - 1)- 
1]. /C 4 (x) = /Ct-i(x) + X^_i(x)(g - 1) = K t -i{x)q = q t 
and JV t (s) = (g* - 1)/(<Z - 1) + ? - 1. 

Therefore the clustering for vertex x after t iterations 

is 



C(x) 



2£ t {x) 



M(x)[M(x)-i] 



(6) 



Notice that A/j (x) = A t is precisely the maximum degree 
at step t. 

We see easily that there is a one-to-one correspondence 
between the degree of a vertex and its clustering. In 
general, we will have that for a vertex v of degree k the 
clustering parameter will be 



C(v) 



2(«-l)(* 



fc(fc-l) 



(7) 



Notice that for q = 2 we have C(v) — 2/k as in Ref. [Tl) . 
Using this result, we can compute now the clustering of 
the graph, 



_ = 2( g + l)(g-l)(A t -f) h g 2(g + !)*"*(, - 1)(A 4 - f ) 



A t (A t -l) 



For t > 7 and for any q > 3 we have 
* - 3g - 1 



A,(A, - 1) 



(8) 
(9) 



Thus the clustering is high, and, similarly to the 7 
coefficient of the power law, is tunable by choosing the 
right value of q: in particular, C ranges from | (in the 
case q = 2) to a limit of 1 when q gets large. 

Diameter. We recall that the diameter of a graph is 
the maximum of all distances between any two vertices of 
the graph. Computing the exact diameter of K(q,t) can 
be done analytically and gives the result shown below. 
Here, we present the main ideas of this analysis. 

First, we note that for any step t, the diameter always 
lies between pair of vertices that have just been created at 
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this step. We will call such vertices outervertices. Then, 
we note that, by construction, no two outervertices cre- 
ated at the same step t can be connected. Hence, if Vt is 
an outervertex created at step t, then vt has been con- 
nected to a q clique composed of vertices created at pair- 
wise different steps t\ < ti < ■ • • < t q and consequently, 
ti < t — q. Now we want to know the distance between 
two outervertices u t and vt- The idea is, starting from ut, 
to "get back" to K(q, 1) by jumps from u t to u t - q , then 
to u t _2q, etc. Hence it takes at most about t/q jumps 
to go from u t to K (q, 1). Then it takes at most as many 
jumps to go from K{q, 1) to v t and we conclude that the 
diameter cannot be bigger than (roughly) 2t/q. 

More precisely, the exact formula for the diameter of 
K(q,t), denoted Diam[K(q, t)], is the following : 

Diam[K(q,t)]=2(\^\+l\+f(q,t) (10) 

where f(q,t) = if t - [(t-2)/q\q < \(q + l)/2], and 
1 otherwise. This value can be obtained by a sharp ana- 
lytical proof. 

When t gets large, then Diam[K(q,t)] ~ 2t/q, while 
Nt ~ g* -1 , thus the diameter clearly grows logarithmi- 
cally with the number of vertices. 

Discussion. The recursive clique-trees K(q,t) which 
we study in this paper are recursively constructed graphs 
which have both small-world and scale-free characteris- 
tics. By choosing adequately the value of q it is possible 



to obtain different clustering parameters and power-law 
exponents (the clustering of the graph ranges from 0.8 
to 1 while the power-law exponent takes values between 
2 and 2.58496). It has been shown that many networks 
modeling complex systems such as the WWW, Internet, 
movie actors, the Erdos number, interactions in protein 
complexes of Saccharomyces cerevisae, the metabolic net- 
work of Escherichia coli, etc., have clustering and power- 
law parameters prec isely in these ranges, see, for exam- 
ple, Refs. 0, H lH [H H3? ]. Moreover, the recursive 
clique trees are actually a deterministic tunable general- 
ization of the scale-free growing networks introduced in 
Ref. |2l| in which one vertex is created per unit time 
and connects to both the ends of a randomly chosen 
edge. These networks were further generalized and called 
"pseudofractal" graph in [11| , and correspond to the par- 
ticular case q = 2 of the recursive clique trees. Variations 
of our construction, for example, by choosing that in dif- 
ferent steps the new vertices added are attached to cliques 
of different size, would allow a richer structure and more 
flexibility in the control of the clustering and power-law 
exponent and other relevant network parameters. 
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